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Bound*  for  Optimal  Confidence 
Limit*  for  Series  Systems 


Abstract 


Lindst ron-Hadden  type  approximations  to  the  lover 
confidence  licit  on  the  reliability  of  a  series  system 
are  theoretically  justified  by  extending  and  simplify¬ 
ing  the  results  of  Sudakov  (1973).  Applications  are 
cade  to  Johns  (1976)  and  Wlnt erbottom  (1974). 


Numerical  examples  are  presented. 
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1.  Introduction  and  Summary 


1  problem  of  fundamental  lnteraet  to  practltlonera  In 
reliability  la  the  statlatlcal  eatlmatlon  of  the  reliability  of  a 
ayates  ualng  experimental  data  collected  on  aubayatcaa*  la  thia 
paper,  the  aubayatem  data  available  conalata  of  a  sequence  of 
Bernoulli  trlale  In  which  a  "one"  la  recorded  If  the  aubayatem 
functlona  ana  a  itro  la  recorded  If  the  aubayatem  falle.  Thue 
for  each  of  "the  k  aubayatema  compoalng  the  ayaten,  the  data  pro¬ 
vided  conalata  of  the  pair  (n^.T^).  1“1. 2, ...k,  where  1^  la 
blnomially  distributed  (n^,p^).  He  aaaume  that  are 

mutually  Independent  random  varlablea. 

The  magnitude  of  lntereat  la  thla  problem  la  eaally  evidenced 
by  the  extenalve  literature  devoted  to  It.  In  thla  regard,  eee 
the  survey  paper  by  Harris  (1977)  and  Section  10.4  of  the  book  by 
Maun,  Schafer,  and  Slngpurwalla  (1974).  In  addition,  the  Defenae 
Advanced  keeearch  Projects  Agency  baa  recently  laaued  a  Handbook 
for  the  Calculation  of  Lower  Statlatlcal  Confidence  Bounda  on 
Syetem  (reliability  (1980). 

Hlatorlcally ,  the  flrat  algniflcant  work  on  thla  problem  waa 
produced  by  Buehler  (1957).  However,  Buehler*e  method  aa  des¬ 
cribed  In  that  paper  la  difficult  to  Implement  computationally 
when  k>2. 

He  proceed  by  describing  Suebler'a  method  In  Section  2.  In 
Section  3  we  apeclallxe  to  series  ayatema,  that  la.  a  ayatem  which 
falle  whenever  at  least  one  aubayatem  fella.  Sudakov'e  (1974) 
reaulta  are  extended  In  Section  4  and  employed  to  exhibit  some 
optimality  propertied  of  the  Llndatrom-Maddcn  method  (aee  Lloyd 
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ar.d  Llpow  (1962))  for  constructing  lower  confidence  bound*  for 
the  reliability  of  serf**  sy*t«a*  of  stochastically  Independent 
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subsystems.  Sone  nunerlcel  examples  ere  given  In  Section  S  end 
the  results  needed  for  this  generalisation  of  Sudakov'a  Theorem 
are  provided  In  the  Appendix  to  this  paper. 

2.  Buebler's  Method  for  Lower  Confidence  Bound* 

A  system  coaposed  of  !c  Independent  eubsyatene  la  said  to  be 
a  coherent  system  (with  respect  to  the  specified  decomposition 
into  subsystems),  if  Che  syscea  fails  when  all  subsystems  fall 
and  che  system  functions  when  all  subsystems  function;  and  replac¬ 
ing  a  defective  subsystec  by  a  functioning  subsystem  can  uot 
cause  a  functioning  system  to  fall.  Coherent  systems  are  des¬ 
cribed  in  Blrabaum,  Esary  and  Saunders  (1961)  and  Barlow  and 
!'i  a. chan  (1975). 

To  any  system  one  can  associate  a  function,  h(p)  - 
‘■\P1»P2‘  •  •  •  »Pg)  •  0<Pl<l.  1-1.2, ...,k.  where  h(p)  Is  the  reliabil¬ 
ity  of  the  system  when  p^  Is  che  probability  that  the  i— 
subsystem  functions.  It  is  well-known  that  If  the  system  Is 

coherent , 

0  <  h(p)  <  1  . 

h(0, • . . ,0)  -  0.  h(l . 1)  -  1  . 

end  h(Pj . Pk)  Is  non-decreasing  la  each  variable. 

For  coherent  systems,  Buehler's  method  may  be  described  as 

follows:  The  observed  outcome  (y,....,y.  )  can  assume  any  of 

k  i  k 

.N  -  II  (n,+l)  values,  since  y.  -  0,1... -,n  . 

i-1  1  1 

vu  denote  n^-y^  by  i-l,2,...,k. 


For  convenience, 


la  said  to  bf  a  noaotonlc  partition,  that  la,  11 

<0,0 . 0)  e  Aj,  (n1,n2,....nfc)  c  A#  and  If  -  (Xjj . »lk>. 

*2  "  <*21.* with  *11  -  *21*  1“1»* . k»  the*  *1  c  *1 

lnpllaa  *2  C  ' AJ *  J  i  *- 

Let 


f(*SP)  “  Pj(X“5> 


k 

a 

1-1 


12.1) 


and  for  1  ^  n  <  fl«  l«t 


lnf|h(p){  l  f(i.ip)  -  o} 
1  jxjCAj.lfn  1  > 


(2.2) 


And  a  •  0. 

S 

Each  such  partition  nay  b«  identified  with  a  function  defined  on 
the  set  oi  sample  outcomes  by  defining  the  ordering  function 
g(x),  where 


g<£)  -  n  if  x  e  Aa  ,  l<n<s  i 


(2.3) 


obviously  g(i)  inherits  the  nonotonlclty  properties  of  the 
partition. 

Subsequently  it  will  be  convenient  to  use  ordering  functions 
g(x)  suct^  that  the  range  of  g(x)  will  be  a  finite  set  of  real 
numbers,  r^<r^< . • • <r# -  With  no  loss  of  generality,  ve  cen  Identify 
the  sets  A^  by  defining  A^  •  {»!*<«)  -  1*1, 2, • . • , s.  We  cen 

new  establish  the  following  theorems. 

Theorem  2.1.  Let  2  be  distributed  by  (2.1).  Then  e  is  e 

l(xj 

(1-c)  lower  confidence  bound  for  h(p).  If  bl(i)  la  also  a  <l-o) 
lower  conf ldtncc  bound  for  b(p),  than  b^<a^,  1<1<s. 


Proof:  Fix  p  and  let  a(p)  ba  tba  a. allcat  Inteaer  such  that 


.fmrrprims «*«i**w*r  > , 


f  n(P>  } 
pp{5  e  A  M 

P-(x  e  u  A,1  >  1-a  . 

Pi  i-n(p)  - 

D»  ‘  fcl’pfc  -  J/i}  *  “} 

-a  confidence  sec  for  p,  sin 

Pp{>  e  Dg(x)}  "  pp{*(i>  i  “<*>}  ^  1-“ 


>  a  . 


and 


Let 


Then  D  is  a  1-a  confidence  set  for  p,  since 

glX) 


(2.*) 


;--s> 


(2.6) 


(2-7) 


Thi-  establishes  the  first  pact  of  the  conclusion.  Further*  since 

h(p)  is  continuous  and  Ofp^l.  the  infleum  in  (2.2)  is  attained. 

How  assume  that  i.  is  the  smallest  index  such  that  b.  >a  . 

x  ix  x1 


X<i  <s-l.  Then,  for  some  5,  p,  , 
-  I-  ox 


b  >  infih(p)  |  l  t  (x;p)  »  a)  -  h(p  )  , 

1  1  XjCAi.ifij  >  ° 

and 


X  f(x;px)  >  a,  h(p1)  <  bt 


Therefore 


p5  (h{5i)  <  b4<i)}  -  -  I  f<*!Pi>  * 
pil  x  t  ixij  «iEA1.i<i1  A 


a  contradiction. 


f t.^e,rk.  Let  d  •  sup-il-h(p))  J  f(x.;p)  -  a>.  Then  4  is  « 
3  1  |x1EAi.i<n  1  >  “ 

(1-a)  upper  confidence  bound  for  l-b(p),  the  unrs 1 lability . 

Let  A  -  |xeEj_.  OtXjCtj,  i-l,2.....k|  and  let  *(x)  be  contin¬ 
uous  on  A  (tbe  closure  of  A)  and  strictly  increasing  in  sscb 


variable  for  x£A.  g(x)  is  to  be  regarded  as  an  ordering  function 


as  described  inncdlately  preceding  Theorem  2.1*  Ve  require  the 
following  additional  property  of  g(x). 

Fix  xq£A.  hat  8(*Q)  <  SU^*0*****^  “  g^-  Then 
g (y ^ • 0, . . . ,0)  -  g(ifl)  has  a  unique  solution  in  y ^ .  Proceeding 


recursively,  let  1^  < 

y,  and 

define  y^  ■  yj(i^)  as  the  solution 

of  *(*„)  “  °* ' 

...,0). 

For  each  l<)fk  and  i^j  <  Fj_j» 

y1>  let 

)  -2  -  -2*  ‘  ’  1  - 

7J  7j '1*  2*  * ' 

solution  of 


g\xo)  ■  g  (  i^  •  ij  •  •  •  •  t  ij  •  y^  •  0  •  •  •  *  »  0)  •  (2  <8) 

We  require  that  the  equations  indicated  la  (2.8)  have 
unique  solutions  for  each  y^ . 

Then  define 

ly2)  Ir2l  lykJ 

fCi  ;p)  -l  l  ...  I  f(isp)  .  (2.9) 

lj-O  i2-o  ifc-0 

where,  for  J>1,  y^  ■  . 

■  sup  F (*  : p)  ,  0<e< J  .  (2.10) 

b(p)-a 

Then  we  have 

Theoren  2.2.  If  5  satisfice  inf  f*(x  ;a)  -  0,  sup  f*(*  ;a)  -  1 

0<a<l  °  0<a<l  ° 

and  f  (*o;a)  is  a  strictly  increasing  function  of  a,  and  if 
xo  e  *n  where  g(x)  detersines  (A2 , *2 , . . . ,A#) ,  and  if 

b  -  inf jh(p) I  l  f(5  ,p)  -  ol  ,  (2.11) 

'  |«eAi,l<n  x  > 

f*(»BJb)  •  a  . 


then  we  have 


Proof;  Since  the  lnflmum  In  (2.11)  In  attained,  there  le  a  Po 

-  ~  ~  *  _ 
such  that  b  -  h(po>  and  F(xo;po>  «  O.  Then  f  (x0,b)  >  a.  It 

r’'(ao>b)  >  a,  there  exist s  p^,  with  a  “  h(?4),  a<b  and 

;>)  -  a  contradicting  (2.11). 

Obviously,  the  above  discussion  can  easily  be  modified  to 
obtain  upper  confidence  bounds  on  the  unreliability  l-b(p)  by 
replacing  inf  by  sup  In  (2.11)  and  requiring  that  f*(xo;a)  be  a 
strictly  decreasing  function  of  a,  0<a<l. 

3.  Applications  to  Series  Systems 
It 

For  a  series  system  h(p)  ■  H  p,.  Further,  throughout  this 

1-1  1 

section  we  assume  that  g(x)  satisfies  the  conditions  necessary  to 
Insure  that  the  solutions  far  y^,...,/-  indicated  In  (2.8)  are 

unique.  Then  we  have  che  following  theorem. 

k 

T.tuorum  3.1.  If  h(p)  -  Dp,,  then  inf  f  (x  ;a)  -  0, 

„  i-1  1  0<e<l  ° 

sup  f  (a  .a)  -  1  and  f*(S  ;a)  Is  strictly  Increasing  la  a„ 

0<a<l  °  ° 

'.hn ever  *  <*ol . *ok>  *OJ  <  W.2 . k. 

Proof .  Since  h(p)  -  1  If  and  only  if  pt  -  1,  1-1,2, ...,k.  It 
follovs  froc  (2.1)  that 

lin  sup  F(x  jp)  -  1  . 
a-»l  h(p)-a 

Similarly,  h(p)  -  0  if  and  only  If  at  least  one  p^  -  0, 

. k.  Since  F(xo;p)  <  Pgjx^nJ  -  1-P-jx^nJ  -  l-q/. 

v <2  have 

llm  sup  1(5  jf)  -  0  . 
a-»0  h(p)«a 

To  show  that  f*(xo;a)  is  strictly  increasing  la  a,  consider 


pak)  satisfy  f  <*„;•)  -  r 


0<a<b<l  and  let  “  (p^ 


Similarly,  let  aatlaly  t  •  F(xo;i>b).  Let 


-  {lx.J 


,  1  >  b«  *ny  non-empty  set  of  indices  such  tb«t 


(~}^^C  <  1  and  let  IC  be  the  remaining  ladlcea.  Then 
t  1/r 

(  a  P  <£>  )  a  p  -  b  . 

jei  *  jel«  alj 

From  the  monotone  likelihood  ratio  property  of  the  binomial 


distribution. 


F(*0JPa)  <  *<*„;?  )  . 


where  the  components  of  p  era  given  by  (3.1).  Then 


F(S0JP  )  *  sup  »<*0!P)  -  P(«o;pb)  -  f  («0!b)  . 

h (p)*b 


4.  Sudakov'a  Method 


f'*>  "  I<r7e7  C  • 


Then  If  y  Is  an  Integer,  >'<a,  ve  have 


!  (?)  p""V  -  I  (n-y.y+1)  . 
1-0  1  P 


For  0<y <u,  real,  define  u(o,y,a)  by  a  -  I  .  . (n-y.y+1). 

Thus,  for  Integer  values  of  y,  u(n.y.a)  Is  e  100(l-a)  percent 

lower  confidence  limit  for  p.  Sudakov  (1973)  showed  that  for 

k 

o,<n  <...<n  and  g(x)  -  fl  (u.-x.i, 

1/  x  1-1  1  1 


•  Ua.p1.«)  <  b  <  u(n1,[y1).a) 


* 


k 

vhcrt  *•  n^qo*  <l0  ” 1 "  i!lit<“i",tot>/*l>  * 

uln^.y^.a)  is  called  the  Lindscroa-tlsdden  method  foe  deteralmlag 
lover  confidence  liilu  for  chc  reliability  of  series  syscess 
(see  Lloyd  end  Lipov  (1962)). 

Lipow  and  Riley  (1959)  used  e  differesc  ordering  function; 
r.  .vorchelcss  they  noced  thsc  for  “snail"  a^,  their  tehuleted 
v.lues  provided  good  sgreeaeat  with  the  results  usiag  the 
Licdstroc-Maddra  method.  For  large  values  of  a^,  the  tehuleted 
values  that  char  provided  are  based  oa  the  Lladstroa-Haddea 
net  hod •  here  we  provide  a  further  Justification  for  the 
Lindstroo-Madden  method  by  establishing  that  it  provides  conserva¬ 
tive  lower  confidence  Halts  (i.e.  is  a  lower  bound  to  b  defined 
In  (2.9))  using  the  ordering  function  |(x)  eaployed  by  Cudekov 
a::d  vr  also  obtain  an  upper  sound  for  b,  thua  deceralnlng  the 
possible  error  of  the  Llndstroa-Maddea  aethod. 

Sudakov's  proof  Is  unnecessarily  complicated  and  contains 
aoae  incorrect  assertions,  which  nevertheless  da  not  affect  the 
validity  of  the  conclusion.  In  the  Appendix  we  provide  a  eiapler 
proof  of  soae  auxiliary  results  needed  for  the  generalisation  of 
Sudakov's  theorem  given  below. 

Theorem  4.1.  Let  g(x)  satisfy  the  hypothesis  of  Theorea  3.1. 

Then, 

b  <  aln  u(a  ,{y*),a)  .  (4.1) 

l<l<k  1  x 

where  b  is  given  by  (2.11)  and  y*  -  3*0*. J2» •  •  •  *Jg_x>  ia  avalnated 
at  J^-0,  2-1,2, . ..,1-1.  Vote  that  y^  -  y*.  If  we  also  have 


y.-l.  yl+l 

>  -iti  .  j-1.2 . k-1  , 

vS "  “j+i 


Proof :  (4.1)  1*  limef  Uu  from  (2.11)  upon  sect  las  fj*l,  J  f  i  sa4 

solving  F(xo;l . 1 . pt , 1 , . . . ,1)  ■  a.  Recall  Chet  »j  <  n2  <  ...  <  nfc 


A. 

?<*o;p)  “  I  b(n1-l1;p1,n1) 


fyk_il 


•  E  b<nL-l"ik-l;pk-l*“k-l>Ipt<"k_,yk**,yk)+1) 

xk-l  u  * 

Now,  apply  L ena«t  Al,  A2,  and  A3  to  the  inatcaott  aua  In  (4.4),  to  got 
Iyk-1* 

£  o  b<“k-l"1k-l!Pk-l,“k-l}rp.  <“k",yk,',yk,+1>  - 
1k-l  * 

lyk-l] 

E  b^“k-l~1k-l:,’k-l,,>k-l^Ip  ^*k"yk,yk+1^  - 

s.-r°  k 

lyk-l> 

Z  b*nk-l~1k-l:,,k-l'"k-l*1p  ^ak«i”Pk-l #yk-l~1k-l+^^  - 

1k-l”°  k 

lpk  1Pk<°k-l"yk-l,yk-l'*'1>  * 

Repeated  applications  of  the  above  sstablisb  Cbat 

P' C*0 » p>  <  i  k  <"l-yl.yl+l)  .  (*.5) 


(4.3)  follows  laaedlatelp  fron  (4.3),  completing  tbs  proof. 
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P.£=orki.  If  (4.3)  bolds  sod  y^  Is  s»  Integer.  tbss  b  - 

It  k<  often  been  suggested  (Lloyd  sad  Llpow  (1M2), 
Kiaterboctoa  (1974),  Bolsbev  sad  Loginov (1966)  ,  Mlraly  sad 
Solo v  *  y lv  (1964))  tbst  tbs  coafidsacs  level  should  depend  ealy  ea 
n,  the  eaallest  sanple  size.  Ke  now  provide  s  numerical  illus¬ 
tration  to  show  tbst  tbe  bound  la  (4.1)  aay  be  laproved  by  cakiag 
all  tbe  a  's  Into  consideration. 

lac  k-3,  a».l,  n  -  (10,12,30),  J  -  (0,3.0).  Thea  far 
3 

S(i)  -  P  -  -341.  «a2.(y2).a)  -  .525. 

t(n^(y^],a)  -  .639.  Tbe  use  of  (4.3)  establishes  .500  <  b  <  .525. 

Mote  tbst  If  *ol  -  n2,  for  sons  1.  l<i<k.  then  |(r)  “  0  sad 
t«0.  It  seems  ressonsble  to  use  b-0  as  the  lower  coafldeace 
Unit  whenever  xoi'n1  for  any  aonotoae  ordering  function  satisfy¬ 
ing  tbe  conditions  of  Section  2. 

k 

Me  now  show  that  If  g(x)  -  I  (n.-x  ),  thea  (4.2)  Is  satls- 

1-1  4  1 

fled  and  Theorem  4.1  applies.  This  result  will  extend  a  result 
due  to  Ulnterboctom  (1974),  who  established  this  fact  for  r 
particular  special  cases.  In  addition,  we  will  also  slow  that 
(4.2)  holds  for  a  number  of  other  ordering  functions  used  in  the 

literature. 

k 

•if.ccrem  4.2.  Let  g(x)  “  D  (n.-x.+a.),  where  «.  >3  »nd 

1-1 

>  a1+1«l(  1-1,2, ... ,  k-1.  Thea  (4.2)  is  satisfied. 

Hroof .  If 

k 

(n.-y.+o,)  n  (n.«».)  -  c 
111  J.tV  J  J 


k 

<“i-k*+«i><“i+i-yiei**i+i)J.;+Jt<aj*n*  - c  • 


ao4 
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then  -t  have 

establishing 

yrki  _  *Wni+ai-ki> 

“i'ki  “i+i  tai+i+“i+i><“i“ki> 


this  last  inequality  will  ba  true  whenever  ni+iai  -  ai+lni"  ln 
particular*  tbls  is  valid  whan  a^-0,  1  »  l*...,k  which  is 
Sudakov'e  ordeclnj  function. 

k 

Theorem  4.3.  If  g(x)  •  1  -  J  x./n.,  than  (4.2)  is  satisfied. 

i-1  1  1 

y  i4-i 

Proof.  If  1  -  y./n,  •  c  *  1  -  - —  -  - *  than 

*  1  “i  *1+1 


or 


This  typa  of  ordering  function  has  been  employed  by  Pavlov 
(1973)*  for  example. 

Theorem  4.4.  Let  g(x)  -  T  «.*.+«  (•?«.)*,  where  s  satisfies 

**o*i  a 

l~$(x  )  *  o  sod  ♦ (x)  !•  the  standard  normal  distribution  function, 
°  k 

•j  >  «2  > - >  «k»  *“*  »j  *  (»t  £  l/®t)  .  Tb<  a  g(i)  satisfies 


(4.2)  if  and  only  if 


li 


<VaJ+l)yJ  -  Ui"*J+l>*a+aJkj"*4kJUa*J+2c'*i<,J*kJ>>  *  <4*4> 


I'roof i  If  g(»0) 


J-l  1=1  * 

c  +  l  tjkj,  then  itfial*!  I  ijkj  * 


i-1 


1-1 


and 


* i7i 


*JkJ 


+  *J+lyJ+l  +  *a(ci+*Jkl+*J+l,)+l> 


*-c  . 


(*.  7) 

(4.6) 


trusting  tb«  left  hand  cldc*  of  (4.7)  and  (4.4),  hi  iktali  (4.4). 
If  k>2,  (4.6)  hold*  for  all  cases  of  interest. 

If  (4.6)  holds,  Chan  setting 


I-a  -  (r(*)>  sj  tx  le  cdt  . 

0 

a  straightforward  Halting  arguaant  shows  that 


wax  a^f  (ly1J-»-l.l-a)  <  h  <  Sjf (jfj+l.l-ai)  .  (4.*) 

Inis  ordering  function  has  bean  uaad  by  Johns  (1976)  and  Vila 
(4.7)  is  the  value  tabulated  by  Johns  for  k«2.  The  validity  of 
the  lower  bound  docs  not  depend  on  (4.6).  In  Table  1  below,  the 
lower  and  upper  bounds  given  In  (4.9)  are  tabulated  along  with 
inc  values  given  by  Johns  for  e».l.  These  refer  to  upper  confi¬ 
dence  Halts  for  the  Poisson  paraacter  combinations  s^A^+SjAj. 

Note  in  particular  that  three  of  the  values  tabulated  by 
Johoc  (Indicated  by  asterisks)  violate  (4.9).  Specifically  consider 
5.24,  la  which  case  (y^)  •  S,  since  g*(2,S)  •  4.7S,  g*(5,0)  •  4.72 
and  g*(6,0)  -  5.46.  Using  the  Poisson  approxlaatlon  we  obtain  the 
value  9.275  for  the  upper  confidence  Halt  to  A  far  a-. I  and  thus 
a^A^+a^Xj  -  5.56.  Consequently  the  sup  suet  exceed  5.56.  An 


alternative  approach  to  the  one  suggested  by  Johns  for  k  >  3  la  to 
alnply  usa  Sjf Oj+l.l-a)  for  b. 

Table  1 


Comparison  of  Upper  ana  Lover  Rounds 
With  Values  Tabulated  by  Johns  for  a*.l 


‘1 

X1 

*2 

Lower 

Bound 

Upper 

Bound 

Johns  * 

Tabled  Value 

.9 

7 

2 

4.79 

5.50 

5.17 

.9 

3 

0 

2.07 

2.27 

2.16 

.75 

6 

3 

6.00 

6.65 

6.23 

.75 

12 

3 

7.90 

8.29 

7.91 

.67 

3 

3 

5.36 

5.61 

5 . 33* 

.67 

15 

2 

8.71 

9.24 

8.81 

.60 

5 

2 

5.56 

5.62 

5.24* 

.60 

7 

6 

9.24 

9.53 

9.18* 

5.  Numerical  Examples  and 

Concluding 

Keaerks 

Examples  1  and  2  Illustrate  tba 

method 

we 

have  described  in 

this  paper. 

k 

Example  1:  Let  H(x)  «  n  (n.-x.),  a 

-  .05. 

k  - 

5.  i  - 

1-1  ^  1 

(20,30,40,25,60),  x  -  (2, 6,10, S, IS) • 

Then 

the 

95X  upper  confl- 

deuce  limit  for  the  failure  probability  la  contained  In  (.86,. 8$). 

Example  2;  Let  H(x)  -  lUn^-x^),  o  -  .05,  k  -  2 ,  5  -  (10,10). 
x  -  (3,2).  Then  the  9SX  upper  confidence  limit  for  the  failure 
probability  is  contained  In  (.70,  .73).  The  value  given  In  Lipow 
and  Riley  (1959)  Is  .70. 

Remarks.  In  this  paper  va  have  shoved  that  the  Llndstrom-Madden 
technique  Is  conservative  for  ordering  functions  satisfying  (4.2). 


Fur t ho 
exact. 
(1974) 
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,  if  y,  is  it  integer,  then  the  Llndstroa-Madden  nethod  is 
Wo  have  also  relaxed  tbs  conditions  needed  in  Wioterbstton 
and  provided  sn  alternative  to  the  netbod  of  Johns  (1*74). 


I 


The  auxiliary  result*  employed  la  the  proof  of  Thtorca  4.1 
arc  provided  here. 

Icama  Al:  I^(n-x,x+l),  0<y<l,  it  a  dccrtttiaf  function  of  n  and 
an  Increasing  function  of  x.  Xy(np#nq+1)#  -  1,  0<p<l,  1 a  an 

increasing  function  of  q. 

Proof :  The  proof  la  laacdiata  from  the  observation  that  the 

beta  distribution  with  peraaetera  a  and  $  baa  aonotone  likelihood 
ratio  in  a  and  -8  and  that  If  a  probability  distribution  baa 
monotone  likelihood  ratl<  in  8#  F^(x)  ia  a  decreasing  function 
of  6  (Lehmann  (1959)*  p.  68  and  p.  74). 


Leprae  A2 :  If 


Xt”^!  ^i+1 

rtir  t  *nd  *<  5  fti+i  •  then 

ni  ki  °i+i  1  1+1 


iy<»1-y1.y1-k1+i)  >  *y<“1+1-y1+1.yt+i+i>  • 

Proof :  ttvrltlot  the  left  and  right  hand  aides  of  (A.l)  aa 

',[<■, -ki><1  -  : 

Leans  Al  appllea  and  the  conclusion  follows. 

Leans  A3:  Let  y.yj  “  y  •  0<y.<l,  1-1,2*  Then 


i*i 

I  (n-x.x+i)  >  l  Kn-kijr  ,n.)I  (n-r.x-k+1)  . 

yly2  k-0  1  1  J2 


U  3) 


Proof: 


1* 


U1  .H.  n-k..  .k  ru-k+l)  P2  _n-x-l„ 

E  <k)yi  (1‘yi>  r(o-x)f<*-k+iT  J  *  ao  4* 

k-o  0 

r(n+l)  U~yX)  yl  fyly2  ft  t yl~t\ X_k  it 

-  T(u-x)  k£0  k!r<x-k+l)  J  l7j7  V  Jx  )  Tj 

r<«*i,  fyxy2  <?> 

-  l-(a-x)  j  kf0  kll'tx-k+lj 


4t  . 


Iku»  (A. 3)  will  bo 14  whenever 


r(nj-X) 

r(n-x>r{x+x) 


y,y2  t""*_1a-t)*4t  > 


p 


rjntil  (yly2  ‘f1  a~yl)^t°~aC'1<yl~t)X"k  gt 
T(a-x)  J  k“0  k!  r(x-k+lj 


or 


rtn-t-l)  fyly2 
lT(n-x)r(x+l)  J 
0 


t»-x-l(1.t)» 


(- 


k-0 


rix±n  /^k/^Cl\x'k]4t  >  0  . 

k:f<x-k+i)  U-t  /  u-t  /  J 


(A. A) 


y.-t  /  l-y,  \ 

Pricing  -  ll  -  ,nd  noting  Chet  O^tty^U.  t<yj  »n4 

{l-t>>l-y1,  we  observe  tbet  (A. 4)  holds  sn4  the  l«w  is  prove4. 
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